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ABSTRACT 

We examine the dependence of the spatial two-point correlation function of quasars 
^qg(r, z) at different redshifts on the initial power spectrum in flat cosmological models. 
Quasars and other elements of the large-scale structure of the universe are supposed to 
form in the peaks of the scalar Gaussian field of density fluctuations of appropriate scales. 
Quasars are considered as a manifestation of short-term active processes at the centers of 
these fluctuations; such processes set in when dark matter counterflows and a shock wave 
appear in the gas. We propose a method for calculating the correlation function iqq{r, z) 
and show its amplitude and slope to depend on the shape of the initial power spectrum 
and the scale R of the fluctuations in which quasars are formed. We demonstrate that 
in the CDM models with the initial power spectrum slope n = 0.7 -h 1 it is possible to 
explain, by choosing appropriate values of R, how the amplitudes and correlation radii 
of Cqq{r, z) may either increase or decrease with increasing redshift z. In particular, 
the correlation radii of ^qq{r, z) grow from 6 — 10h~^ Mpc when R grows from 0A5h~^ 
to 1.3/i~^ Mpc. The H+CDM model at realistic values of R fails to account for the 
observational data according to which the (^qq(r, z) amphtude decreases with increasing 
z. 

Subject headings: cosmology: initial power spectrum - dark matter - quasars 



1. Introduction 

The most elaborated scenario of the origin of the 
large-scale structure in the universe is pictured by the 
model in which such a structure results from the evo- 
lution of the uniform isotropic Gaussian scalar field 
of cosmological matter density fluctuations under the 
effect of gravitational instability. The major unre- 
solved problems in this case are the choice of an in- 
flationary model for an early universe, determination 
of the nature of the "dark" non-baryon matter (DM) 
and its fraction in the total matter mass. Assump- 
tions as to the inflationary model and DM nature 
define the shape of the initial (post-recombination) 
power spectrum of cosmological density fluctuations, 
P(fc), and thus the principal parameters of large- 
scale structure can be theoretically calculated. There- 
fore it is of importance to test cosmological mod- 
els with given power spectra P{k). This testing can 
be done by calculating spatial two-point correlation 
functions for large-scale structure elements on var- 
ious scales and comparing them with observational 
data. The testing is based on the relationship be- 
tween the characteristics of structure elements and 
their correlation functions, on the one hand, and 
the amplitude and slope of the initial spectrum P{k) 
at different fc, on the other hand. Information on 
the power spectrum on small and intermediate scales 
(r < lOO/i-l Mpc, h = Ho/ilOOkm ■ s'^Mpc-^, Hq 
being the Hubble constant at the present epoch) re- 
sides in the correlation functions of bright massive 
galaxies, rich clusters of galaxies, and quasars. The 
"observed" correlation functions of all these three 
types of objects are described by the approximate 
expression ^(r) = (r/rp)"^'®, where tq is the corre- 
lation radius equal to 5Ah-^ Mpc, 18h-^ Mpc 
(6-=- 10)h-^ Mpc ||-||. for the three types of objects, 
respectively. At the moment investigations of the 
" observed" correlation function of quasars in different 
redshift ranges give ambiguous results. For example, 
the correlation function amplitude found in ^ |^ 
decreases with increasing z, it remains unchanged in 
1 4], and grows at z < 1.7 and diminishes at z > 1.7 in 
7]. The common result of these studies is that the 
amplitude for quasars is larger than for galaxies but is 
smaller than for rich clusters of galaxies. This result 
was obtained by different authors from two quasar 
samples: a combined sample of all quasars observed in 
the z range from 0.1 to 4.5 and the "nearby" quasars 
at z < 1.5. 



Here we look into the possibility of using the above 
results as a test for cosmological models with given 
power spectra P{k); to this end, the theoretical cor- 
relation function of quasars has to be calculated. 

Theoretical methods for calculating correlation func- 
tions of galaxies and their clusters are based on the 
theory of Gaussian random fields [p[-|ll|], they have 
been devised in detail. The results obtained within 
the scope of cosmological models with given initial 
power spectra were analyzed in Jl2[]l^ , for example. 
Calculation of the correlation function of quasars 
is complicated by a number of problems which have 
yet to be resolved. What scale is typical of the re- 
gions where quasars formed at different zl What are 
typical physical parameters of quasars: mass, dura- 
tion of formation, lifetime, etc.? What is the re- 
lation between ^qq and these parameters? Answers 
to these questions essentially depend on the physical 
model chosen for the quasar phenomenon. In par- 
ticular, the disk accretion of gas on a massive black 
hole at the center of a galaxy may be a model mech- 
anism. Therefore, for the mass of the fluctuations in 
which quasars are formed we can take the mass of 
"parent" galaxy Mg/q which is able to ensure a high 
luminosity of the nucleus (observed as a quasar) dur- 
ing the quasar lifetime Tq. Based on the results of 
||]l|, we may take 10" - lO^^M© for Mg/q (the 
black hole mass - 10^ - lO^M©), - 10^ - 10^ yr for 
Tq, and Atq < 10^ yr for the duration of quasar for- 
mation. Whether these parameters are the same for 
quasars at different z is still an open question. 

2. Principal assumptions and formulation of 
the problem 

In the cosmological scenario used by us, galaxies, 
rich clusters of galaxies, and quasars appear in the 
peaks of the scalar Gaussian field of matter density 
fluctuations on corresponding scales, the relative am- 
plitude of fluctuations being 5 = 5p/ p — a ■ v [a is the 
rms amplitude, v is the peak height). It is assumed 
that galaxies and their clusters come into being when 
counterflows arise in the DM and a shock wave arises 
in the gas. The amplitude 5 at this moment t is de- 
termined from Tolmen's model in terms of redshift: 
5{z) = 1.69 • {z + 1). The amplitude corresponding to 
the objects that appeared earlier is 5 > 5{z), it has a 
normal distribution: p{5) = (27rcr^)~^/^ • exp{—5'^/2a'^) 
The probability that a galaxy or a rich cluster of 
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galaxies occurs at a fixed z is 



Expression (6) is simplified at ^(r) <C 1 and \ 



Piiz)^ J piS)dS. 

5{z) 



(1) 



The probability that two galaxies exist simultane- 
ously at a fixed z at two different points xi and X2 
(r = \xi - X2\) is 



P2iz) ^ J J P{Sl,S2)dSidS2, (2) 

5{z)S{z) 

where p{di,S2) is the two-dimensional normal distri- 
bution of random amplitudes 6i and S2 : 

p{d^,62) - (2 • 7:)-' ■ [y/eio,z)-e{r,z)y' 



xexp 



aO,z)-Sl+mz)-Sj-2-^{r,z)-S,-S2 



(3) 

here r > and ^(r, z) is the correlation function of the 
density fiuctuations in which the objects are formed. 
The function is calculated from the given initial power 
spectrum P{k,Rf) smoothed on the scale Rf which 
corresponds to the scale of the objects: 



- I ,2 P{k,Rf) sin{kr) 



1 

2^ 



(i + zy 



kr 



dk, (4) 



where 



and 



P{k,Rf) = P{k) ■ W^ikRf) 



W{kRf) = exp i^-- -k^ -Rj 

is the smoothing function. The statistical correlation 
function of the fluctuation peaks in which cosmologi- 
cal objects are formed is, by definition, 



..t,. .^ _ P2{z) 
P^{z) 



C{r,z) 



1. 



(5) 



This function for rich clusters of galaxies or for galax- 
ies at z = is 



Clir) = U{r,z = 0) = \l--{ erfc ( ^ ) I x 



V 



f e-'/^ y" -erfcl 

1/ ^ 



-e{r)/em, 



dy-l. 
(6) 



(7) 



In this case Coo('') ^'^^ objects is related to ^(r) for 
the fluctuations on the corresponding scale in which 
the objects are formed. The factor v/a \s called the 
statistical biasing of these objects, and the procedures 
for its determination for galaxies and rich clusters of 
galaxies are described in |o|. In dor iving the 

correlation function of quasars we take advantage of 
the approach proposed in [pd| with allowance made 
for the specific nature of quasars. 

3. Correlation function of quasars 

To calculate correlation functions of galaxies and 
rich clusters of galaxies, we have to specify the scale 
of the corresponding fluctuations in which they are 
formed. For quasars, there are two more impor- 
tant quantities along with this scale: the duration 
of quasar formation Atg and quasar lifetime Tq. The 
smallest mass of the galaxies which may experience 
the quasar phase in the course of their evolution is 
^s/g" - 2 • IO^Mq |ll-|i. The duration At, < 10^ 
yr is much shorter than the cosmological evolution 
time of fluctuations in which such galaxies are formed 
at z < 5, and so we may ignore this quantity in first 
approximation. The same studies |T5| - |l8| reveal that 
the quasar phase is short-lived Tq ^ 10^ — 10® yr. The 
quantities Pi{z) and ^2(2) in (1), (2) are the proba- 
bilities that the random amplitudes of corresponding 
fluctuations are 5 > 5{z)^ i. e., they determine the 
probability of coming into being or the probability of 
the existence of objects on corresponding scales with 
redshifts larger than a given z. For quasars, the in- 
terval l^{z,Tq) between the limits of integration in (1) 
and (2) is determined as the lifetime Tq of the quasars 
which came into being at the cosmological moment 
t{z) corresponding to a given z |2ll-p3[: 



When 



A(z,T,) = 1.13- (2 + 1)-^. (8) 
t(z), the quantities Pi{z) and P2{z) are 



in form (1), (2) with the upper limit of integration 
5{z) -f A(z,Tg). When Tq <C t{z), expressions (1), (2) 
take the form 

Pl{z)^p{6{z)y^{z,Tq), Pl{z)^p{5,,52)-^\z,Tq). 

(9) 
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in view of smallness of A(z, r^). Then, on the basis 
of (5) and in view of (9), the correlation function of 
quasars is 



X exp 



2.86- (z + 1) 



(M\' 
\m) 

, i(0) 



1, (10) 



As it follows from (4), the ratio ^(r, z)/^(0, z) is inde- 
pendent of z, and ^(r) = ^(r, z = 0). 

The above expression for ^q*q{r,z) represents the 
statistical component in the correlation function of 
the fluctuation peaks where quasars are formed, but 
it takes no account of the dynamics of background 
large-scale inhomogeneities (the dynamical compo- 
nent). According to the complete correlation 
function of objects is 

eoo(r, z) = (y^^^) + ^C'„(r, z)^ , (11) 

where the statistical component for quasars (,l*o{r, z) = 
^**(r, z) is represented by (10) and the dynamical 
component (^^^(r, z) = ^(r, z) is correlation function 
(4) of density fluctuations. 

Let us examine the approximation (C(^)/C(0)) ^ 1- 
In this approximation £,^„l,{r, z) (10) takes the form 



Cir,z) 



1.69- (1 + z) 



(12) 

after its expansion into the Taylor series. It coincides 
with expression (7) for the correlation function of 
the galaxies and clusters which formed in high peaks 
u ^ 1 (that is, at high redshifts), but in the case of 
quasars there is no such a restriction to high peaks 
only. 

According to (4) , the dynamical component ^ (r, z) 
is expressed in terms of ^(r) 

e(r,z)=f(r).(l + z)-^ (13) 

Then, in view of (11)-(13), we can write 

■1.69- (l-fz) 1 



1 + z 



iir). (14) 



As seen from (12) and (13), the statistical component 
in the correlation function of quasars increases with 



z, while the dynamical component diminishes. There- 
fore, when z is variable and r is fixed, the correlation 
function amplitude may have a minimum if a decrease 
in the dynamical component is not compensated by 
an increase in the statistical component at small z. 
The minimum point can be found from the condition 
d^qq{r,z)/dz — with expression (14): 



a 



1. 



(15) 



Thus, at (7 < 1.3 we have d£,qq{r, z) jdz > within the 
interval z > the correlation function amplitude does 
nothing but grows. When cr > 1.3, the amplitude is 
minimum at z — z^. Over the entire range z < 5, 
in which quasars are observed, the amplitude of the 
correlation function of quasars diminishes at cr > 7.8, 
over the range z < 4 it diminishes at ct > 6.5, and over 
the range z < 3 at cr > 5.2. These results were derived 
in the approximation {£,{r) / S,{0)) ^ 1, which is valid 
at r > 10h~^ Mpc. On these scales, the large-scale 
inhomogeneities have the amplitude < {Sp)/p ><C 1 
even at z 0, so that the correction of the spectrum 
P{k) for the nonlinear evolution does not alfect the 
results. 

4. Calculation results 

We calculated the correlation function of quasars 

5, qq{r, z) within the framework of models with various 
initial power spectra of fluctuations P{k). Among 
models with various (fraction of baryon density in 
the density of the whole matter in terms of the crit- 
ical density fib = Pb/pcr), ^cdm, and Qhdm (CDM 
stands for the cold dark matter of axion type and 
HDM for the hot dark matter of massive neutrino 
type), the CDM models with fli, < 0.1 and flcDM = 
1 — fib are believed to be the most promising at the 
moment; they include the "standard" CDM model 
with the spectrum slope n = 1 (-P(fc) = A - ■T^{k), 
where A is a spectrum normalization constant, k is 
the wave number, and T(fc) is a transfer function de- 
pending on DM nature), "inclined" CDM models with 
n = 0.8; 0.7, and the "hybrid" H-hCDM model with 
ilb = 0.1, flcDM = 0.6, flHDM = 0.3, n = 1. Model 
parameters and methods of spectrum normalization 
on the COBE results (24[ H were described in ^ 
^,|2~ 



To calculate S,qq{r, z), one has to specify the 
characteristic scale or mass of fluctuations in which 
quasars are formed and smooth the power spectrum 
P{k) by a Gaussian filter with the radius Rf corre- 
sponding to this scale. It follows from the expression 
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M = 4.37- W'^Rjh-^MQ |ig that the radius corre- 
sponding to the mass M^™ ~ 2 • lO^Mg IHlS 
Rf ~ 0.35/i~^ Mpc. So, we start from the assump- 
tion that quasars are an early short phase in the evo- 
lution of galaxies with this mass and calculate the 
correlation function of quasars at this value of Rf. 

Figure 1 shows approximations of the observed cor- 
relation functions of galaxies |§, ||] and rich clusters of 
galaxies [|l|, D and the correlation function of quasars 
calculated by expression (11) for z = 0.5,1,2,3,4 in 
various cosmological models. 

It is apparent that at large redshifts, where the 
correlation function is mainly specified by its statis- 
tical component and the dynamical component is of 
minor importance, the amplitudes of correlation func- 
tions of quasars diminish with increasing n in var- 
ious CDM models. This is a manifestation of the 
well-known tendency of Gaussian field peaks to de- 
crease the degree of their clustering with decreasing 
height I/. Calculations of these heights for a fixed z 
by the expression v = 6{z)/a = 1.69 • (1 -t- z)/a where 
a = a{Rf) — ^/^(0) (4), reveal that they are the 
largest in the CDM model with n ~ 0.7 (rr = 1-67), 
somewhat smaller in the CDM model with n = 0.8 
(a = 2.37), and the smallest in the CDM model 
with n = 1 {a = 4.78). The above rms amplitudes 
a allow one to follow the variations of correlation 
function amplitudes for quasars in these models at 
different z (see expression (15)). Thus, the ampli- 
tude ^gq{r, z) in the CDM model with n — 0.7 di- 
minishes at small redshifts, reaches its minimum at 
z z 0.3, and then it grows. In the CDM model 
with n = 0.8, a similar minimum occurs at z « 0.8, 
in the CDM model n ~ 1 at z ^ 2.7, and in the 
H-hCDM (cr = 1.39) - at z « 0.07. These resuhs are 
in good agreement with numerical calculations (Fig. 
1). When the plots in Fig. 1 are matched to ob- 
servational data , one can see that only the CDM 
model with n = 1 and results of |^ are in accord as to 
the correlation radius ro{z = 0.5) ~ (6 10)h~^ Mpc 
and the trend of the ^gg(r, z) amplitude variations. 
As regards other observational data, the cosmologi- 
cal models discussed here under the assumption that 
quasars are formed in the fluctuations on the scale 
Rf = 0.35h~^ Mpc fail to explain the amplitudes of 
the observed correlation functions of quasars at var- 
ious z and their correlation radii. According to 
|, the radii ro(z) ~ (6 10) /i"! Mpc at different red- 
shifts. These correlation function parameters in the 
CDM models are smaller than in the observed func- 
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Fig. 1. — Correlation functions of quasars ^qq{r, z) 
calculated in different models at Rf = Q.35h~^ Mpc 
and correlation functions of galaxies (circles) and rich 
clusters of galaxies (squares) approximated by the ex- 
pression ^(r) = {r/ro)~^-^ at Tq « 5.4/i~^ Mpc and 
rfj' « 18/1-1 Mpc. 
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tion, while in the H+CDM model they are too large. 
This may mean that the fluctuation scale chosen by 
us, Rf = Q.35h~^ Mpc, is too small in the CDM mod- 
els and too large in the H+CDM model. The corre- 
lation functions ^qq{r,z) plotted in Fig. 2 were cal- 
culated with the CDM models with such fluctuation 
scales Rf for each model that the functions do not 
contradict the data of |Q, || ^ on the amplitude (or 
correlation radius) on the scales r < 12h~^ Mpc. At 
12 < r < 40 Mpc the amplitudes of these correla- 
tion functions of quasars are no smaller than the am- 
plitude of the observed correlation function of galax- 
ies, with due regard for a large scatter of that func- 
tion on these scales. Besides, to make more correct 
the comparison between the theoretical and obser- 
vational data, we took into consideration that the 
observed correlation functions of quasars were deter- 
mined for wide redshift intervals Az = [zmin , Zmax] 
1^^, ^, while the theoretical functions were calcu- 
lated for fixed z. Therefore, the observed correlation 
functions at different Az were fitted by the theoret- 
ical functions calculated for z, the weighted mean in 
these intervals found from the observed distribution 
of quasars nq{z) in the comoving reference frame: 

J z ■ n(z) ■ r^(z) dr(z) 

/ n(z) • r^(z) dr[z) 

where r{z) = 2c/ Hq • (l — (z + 1)^" ''') is the distance 
in the comoving reference frame, dr{z) = c/ Hq ■ (z ~\- 
dz, and c is the velocity of light. The weighting 
made to fit £,qq{r, z) to the results of [|[ ^ gave z = 
1.3 for the interval Az = [0.1, 1.5] and z = 2.6 for 
[1.5,4.5]. For the intervals Az used in studies [|[ ^ 
we found z = 1.05 for [0.1, 1.1], z = 1.49 for [1.1, 1.7], 
and z = 2.38 for [1.7,3.1]. 

Based on the correlation functions of quasars which 
correspond to the data of [Q, ^, ^ as to correlation 
radius rQ(2;) ~ 6 lO/i^^ Mpc, we determined the 
scales Rf oi fluctuations in which quasars are formed 
and their rms amplitudes cr{Rf) in the CDM models: 
Rf « 1.3/1-1 Mpc, cr(i?/) = 2.24 (in the model with 
n = 1); Rf « 0.7/1-1 ^p^,^ ^ ^ q g). 

Rf « 0.45/1-1 Mpc, a{Rf) = 1.52 {n = 0.7). 

Judging from the above values of cr{Rf) and ex- 
pression (15), the amplitude and the correlation ra- 
dius of £,qq{r, z) in these models grow with z. 
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Fig. 2. — Correlation functions of quasars 
£^qq{r, z) calculated in CDM models for z = 
1.3, 2.6, 1.05, 1.49, 2.38 weighted means in the Az in- 
tervals 0.1- 1.5, 1.5-4.5, 0.1-1.1, 1.1-1.7, 1.7-3.1. 
The scale Rf is chosen such that the correlation radii 
may be within the range 6 — 10h~^ Mpc. Other des- 
ignations are the same as in Fig. 1. 
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5. Conclusion 

We have elaborated a method for the theoretical 
calculation of the correlation function of quasars, it is 
based on the Gaussian statistics of the initial field of 
cosmological matter density fluctuations. The corre- 
lation functions of quasars were calculated for differ- 
ent rcdshifts in cosmological models with given initial 
power spectra P{k) on the assumption that quasars 
are formed in the peaks of the fluctuations and ex- 
ist over a time much shorter than the cosmological 
time. The amplitudes and slopes of the correlation 
functions are shown to depend on fluctuation scale 
and power P{k) on corresponding scales. The ampli- 
tudes and correlation radii of the observed and theo- 
retical correlation functions of quasars are consistent 
in the CDM models with the spectrum slope n ranging 
from 0.7 to 1 when the scale Rf of the corresponding 
fluctuations in which quasars are formed ranges from 

0. 4.5h~^ Mpc to 1.3h~^ Mpc. In this case, however, 
the function amplitudes increase with z. Although 
there is no evidence of any monotonic growth of the 
correlation function amplitude in studies ^, ^, |^, H] , 
the estimates are so contradictory that we cannot use 
them to test the cosmological models. The H-t-CDM 
model fails to explain the correlation radius derived 
from observational data (6/i^^ ~ lOh^^ Mpc) if the 
scale of fluctuations in which quasars are formed ex- 
ceeds 0.35/i~^ Mpc. Our results are in accord with 
the physical models which regard the quasar phe- 
nomenon as an early short phase in the evolution of 
massive galaxies or as a merger of galaxies situated in 
groups. 
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